Graphene valley filter using a line defect 
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With its two degenerate valleys at the Fermi level, the band structure of graphene provides the 
opportunity to develop unconventional electronic applications. Herein, we show that electron and 
hole quasiparticles in graphene can be filtered according to which valley they occupy without the 
need to introduce confinement. The proposed valley filter is based on scattering off a recently 
observed line defect in graphene. Quantum transport calculations show that the line defect is 
semitransparent and that quasiparticles arriving at the line defect with a high angle of incidence 
are transmitted with a valley polarization near iOO%. 

PACS numbers: 73.22.Pr, 73.61.Wp, 73.63.Bd, 85.75.-d 
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Owing to its exceptional electronjj and thermal[2,] 
transport properties, graphene Q is a promising material 
for use in advanced low-energy electronic applications. 
As graphene is a semimetal with no band gap,[3l it can- 
not in its intrinsic form be used as a replacement mate- 
rial for silicon and other semiconductors in conventional 
electronics. However, with its bands at the Fermi level 
defining two conical valleys, graphene might instead of- 
fer novel electronic applications. Before such applications 
can be realized, the control over transport in graphene 
needs to be improved. Understanding the effects a re- 
cently observed line defect Q and grain boundaries @ 
have on the transport properties might be the key. While 
there are several theoretical studies that have investi- 
gated the latter grain boundaries, the line defect 
has so far remained relatively unexplored, even though 
the line defect is always straight and its adjoining grains 
are aligned, arguably making the line defect a more suit- 
able structure for controlled transport in graphene. 

This Letter investigates the transport properties of the 
atomically precise, self- assembled graphene line defect 
shown in Fig. [I] and shows that this line defect is semi- 
transparent and can be used as a valley filter. Electron 
and hole quasiparticles can either transmit through the 
line defect without changing direction or reflect follow- 
ing the law of specular reflection. The transmission and 
reflection probabilities depend on the valley degree of 
freedom, thus allowing the quasiparticles to be filtered 
according to their valley degree of freedom with a po- 
larization near 100% for quasiparticles arriving at the 
line defect with a high angle of incidence. This fil- 
ter is much different from the valley filter proposed by 
Rycerz et al.,fn\ which relies on the isolation of a few 
one-dimensional channels in a narrowly confined region. 
The latter filter has not yet been demonstrated, presum- 
ably due to challenges in fabricating the structure, which 
requires a sub-lOnm constriction with saturated zigzag 
edges. In contrast, the filter proposed herein relies on the 
two-dimensional geometry of graphene and its required 
structure has already been observed, j5j as demonstrated 
by the micrograph in Fig.[TJ The valley filter is expected 
to be a central component in vallevtronics.il l| just as the 
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FIG. 1. Extended line defect in graphene. (a) Arrangement of 
the carbon atoms around the line defect highlighted in gray. 
The structure exhibits translational symmetry along the de- 
fect with a primitive cell shown in beige. The two sublattices 
in graphene indicated by blue and green atoms reverse upon 
reflection at the line defect, as can be clearly seen in the en- 
larged overlay, (b) Scanning tunneling microscope image of 
the line defect in graphene (Image adapted from Ref. [^). 



spin filter is central in spintronics. Electronics that make 
use of the two valleys in graphene is attractive because 
the valleys are separated by a large wave vector, making 
valley information robust against scattering from slowly 
varying potentials. [l^ including scattering caused by in- 
travalley acoustic phonons that often limit coherent low- 
bias devices to low-temperature operation. Therefore, 
the valley information generated by the filter proposed 
herein could in principle be preserved even in a diffusive 
charge transport regime. 

Consider a low-energy electron (hole) quasiparticle 
with energy s and valley index r approaching the line 
defect from the left (right) at the angle of incidence a. 
Asymptotically far from the line defect, the quasiparticle 
occupies a graphene state |$t), where r = ±1 is a val- 
ley index. See Fig. [51 Let the quasiparticle wave vector 
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FIG. 2. Valley states scattering off the line defect, (a) Ener- 
gies close to the Fermi level in graphene have dark contours 
and are located in the corners of the first Brillouin zone con- 
tained within the gold hexagon. The two valleys, K and K', 
are identified as the two disjoint low-energy regions in the 
reciprocal primitive cell enclosed by the blue rectangle, (b) 
An incident quasiparticle state is defined by the valley index 
r and wave vector g*, where the latter points in the direction 
q given by the angle of incidence a. (c) Owing to energy 
and momentum conservation along the line defect, there are 
only two nonevanescent scattered states allowed, (d) The sub- 
lattice symmetric |-|-) and antisymmetric | — ) components of 
the incident state \'^t) are transmitted and reflected, respec- 
tively. The thickness of each arrow indicates the probability 
the quasiparticle will follow the respective path. 



q — ((/a;, qy) be the wave vector measured from the center 
of the occupied valley, located at Kr ~ AnTy/ia, where 
a is the graphene lattice constant and y is the unit vector 
along the line defect. To first order in g = the nearest- 
neighbor tight-binding Hamiltonian in graphene [3| can 
be expressed [l3j| as Hr = fivp [q-xCy + rqyax), where 
vf = ^/3\^\a/2h is the Fermi velocity with the nearest- 
neighbor hopping parameter 7 w — 2.6cV, and and ay 
are Pauli matrices. The Hamiltonian has energy eigen- 
values E = rje, where e = hvpq and 77 — +1 (—1) if 
the quasiparticle is an electron (hole). From the quasi- 
particle energy dispersion e, it follows that the quasi- 
particle group velocity is vpq, where q is the unit vector 



in the direction q. Because the quasiparticle travels in the 
direction g, q^ = rjq cos a and qy = ijq sin a. Using these 
relations, the eigenstate of the graphene Hamiltonian for 
a given r and a can be expressed as 



(1) 



where \A) and \B) refer to the two sublattices in 
graphene. 

The structure in Fig. [T] exhibits a useful symmetry line 
through the line defect. In the limit g — >■ 0, the reflec- 
tion operator commutes with the graphene translation 
operator perpendicular to the line defect. Therefore, 
symmetry-adapted states |±) can be constructed that 
are simultaneous eigenstates of the graphene Hamilto- 
nian and the reflection operator. As the reflection oper- 
ator maps A sites onto B sites, and vice versa, it can be 
represented by the operator ax acting on the two sublat- 
tices. From the eigenstates of a^, one obtains 
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The graphene state ^ expressed in the symmetry- 
adapted basis is 
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The full Hamiltonian describing the system can be 
divided into three terms, Hr — Ht + H]j + V, which 
represent graphene, the isolated line defect, and the in- 
teraction between graphene and the line defect, respec- 
tively. As each term commutes with the reflection op- 
erator, the full Hamiltonian must commute with the re- 
flection operator, and thus, the eigenstates of Ht in the 
symmetry-adapted basis are either symmetric or anti- 
symmetric about the line defect. Antisymmetric states 
have a node at the line defect, and as a result, there are 
no matrix elements within the nearest-neighbor model 
coupling the left and right sides. Therefore, antisymmet- 
ric states cannot contribute to any transmission across 
the line defect. As shown below, however, there are two 
symmetric states at the Fermi level without a node on 
the line defect. As these states are extended eigenstates 
of the full Hamiltonian, they carry quasiparticles across 
the line defect without scattering. Thus, we can conclude 
that the transmission probability of the quasiparticle ap- 
proaching the line defect is 

Tr = \{+\^r)\^ = I i^ + r sin a). (4) 

As the sum of the transmission probabilities in Eq. (|4]) 
over the two valleys is T,- = 1, we can also conclude 
that the line defect is semitransparent. The semitrans- 
parency follows from the relation (-|-|$_^) = {—\^r)* and 
the normalization of 1$,-). See Fig. [51 Fig. [3] shows that 
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FIG. 3. The probability that an incident quasiparticle at the 
Fermi level with valley index r and angle of incidence a will 
transmit through the line defect. 

the transmission probability of a quasiparticle varies sig- 
nificantly with its angle of incidence a. At a high angle of 
incidence, there is almost full transmission or reflection, 
depending on the valley index r. 

Owing to the semitransparency, given an unpolarized 
beam of incident quasiparticles, the probability Pr that 
a transmitted quasiparticle has valley index r is given 
by Tr. The polarization "P = (t) P+i - P_i of the 
transmitted beam is 

V = sin a. (5) 

This expression shows that an unpolarized beam of quasi- 
particles approaching the line defect at a high angle of 
incidence will lead to outgoing transmitted and reflected 
beams that are almost completely polarized. 

The results presented above are based on symmetry 
arguments that are valid only in the limit as the quasi- 
particle energy e — >■ 0. The results, however, hold to an 
excellent approximation as long as e ^ fivp/a ~ 2.3 eV. 
To show this, we have performed numerical transport 
calculations that are treated exactly within the nearest- 
neighbor tight-binding model. The transmission proba- 
bility is shown in Fig.|3]as a function of the components of 
the wave vector q. Tracing the transmission probability 
in the figure along any constant energy contour, which is 
almost perfectly circular due to the approximately conic 
dispersion, yields a dependence on the angle of incidence 
that is virtually indistinguishable from that in Fig. [31 thus 
confirming that almost full polarization can be achieved 
as long a,s qa <^ 1 ov e <^ hvp / a- To further test the 
robustness of the valley filter, we have performed calcu- 
lations with interactions across the line defect and with 
potentials on the line defect and their neighboring sites. 
None of these tests led to results qualitatively different 
from those presented in Fig. 21 

Because the structure in Fig.[lja) exhibits a symmetry 




[7i/a] 

FIG. 4. Transmission probability for a quasiparticle with fi- 
nite energy, (a) The transmission probability for a quasi- 
particle approaching the line defect at an angle of incidence 
a is indicated by the brightness, ranging from to 1 in in- 
crements of 0.05. For the angle a shown, the transmission is 
near one if the quasiparticle has valley index r = -1-1, as illus- 
trated on the right, where the filled (open) circle represents 
an electron (hole) quasiparticle. (b) The transmission prob- 
ability for the corresponding quasiparticle with valley index 
r = -1. 

plane through the center of the primitive cell, there is 
time-reversal symmetry in the direction along the line 
defect. This time- reversal symmetry implies that the 
transmission probability of a quasiparticle with valley 
index — r can be obtained from that of a quasiparticle 
with valley index r by letting qy — ;> ~<lv This rela- 
tionship between r = ±1 can be seen in Fig. SI Note, 
however, that the transmission probabilities in Fig. SI are 
not symmetric about qy, and thus one can conclude that 
the scattering of a quasiparticle depends on the valley 
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index, which is a necessary requirement for a valley fil- 
ter. As both graphene and the line defect have Hamilto- 
nians exhibiting electron-hole symmetry, one might ex- 
pect that Tr{q) = Tr{—q). The scattering does, in gen- 
eral however, not obey electron-hole symmetry. That 
the condition is satisfied by Eq. is in part a conse- 
quence of the neglect of evanescent waves. These evanes- 
cent waves are accounted for in the numerically obtained 
transmission probability in Fig.Hl Note, for instance, 
that r^(0.0057r/a,0.0l7r/a) 9^ r^(-0.0057r/a, -0.0l7r/a). 
To understand the lack of electron-hole symmetry in the 
combined graphene-line defect system, we note that the 
structure in Fig.[TIa) is not bipartite; in particular, the 
sites participating in the pentagons at the line defect can- 
not be divided into two types where one type has only 
nearest neighbors of the other type. 

To gain further insight into the numerical calculations 
and how they lead to Eq. (HJ, it is useful to perform the 
transport calculations analytically in the limit of small q. 
As there is translational symmetry along the line defect, 
one can focus on those atoms within the primitive cell 
shown in Fig.[T] The Hamiltonian of the isolated line 
defect is then 

Next, we seek a retarded self energy S that accounts 
for the coupling of the two line defect atoms to the semi- 
infinite portion of graphene on each side of the line defect. 
As shown in Fig.[Tl there are two atoms neighboring the 
line defect on each side. Expressed in the basis of the 
atoms parallel to the line defect, the graphene state with 
valley index r is given by 

1^) = g-2ir/3 ^ ■ (7) 

As the basis contains two atoms belonging to the same 
sublattice, the graphene state above can be folded onto 
another graphene state with the same wave vector in the 
full system. This latter state, however, is evanescent near 
the Fermi level and can be neglected. Requiring E to be 
retarded fixes the relative phase between the atoms on 
the line defect and their neighbors, resulting in the rela- 
tion S(i?|$T-) = 7(A|$T-)|r)(T|, from which one obtains 

Equipped with Hd describing the interactions within 
the line defect and E describing the coupling to the 
semi-infinite graphene on each side, one can calculate 
the retarded Green function on the line defect, G — 
(r]el — Hd — 2E) ^, where / is the unit matrix. To ze- 
roth order in q, 

1 -t- ie»^" ^ 1 - ie*^^"-^'^/^) jg^ra j • 

(9) 



The probability that the quasiparticle will transmit 
through the hne defect is given by = {t\VGVG^\t), 
where F = i(E-E^). Inserting Eqs. (0111) into this 
equation, one recovers Eq. exactly. 

In the initial analysis leading to Eq. Q , an assertion 
was made that there are two symmetric states at the 
Fermi level without a node on the line defect. That 
claim can now be verified using Eqs. ([TKnilHl). Accord- 
ing to Eq. ([1}, a symmetric state must satisfy tol — 7r/2. 
When this condition is satisfied, one finds that the deter- 
minant det [Hd -\- 2E) = 0, which implies that there are 
exactly two symmetric states at the Fermi level, one for 
each valley index t. The corresponding eigenstates are 
I^E't-) = I — t), confirming that the states have no node at 
the line defect. 

Rather than forming isolated Bloch waves, experi- 
ments exploiting the valley filter will likely construct 
more complex wave patterns. The dimensions of the 
system should be chosen such that the mean-free path 
is longer than the distance between the source and the 
line defect. For low-energy quasiparticles, the wavelength 
A = 27r/g is much greater than the repeating length 
2o of the line defect. As long as this repeating length 
is also much shorter than any spatial features of the 
waves, to an excellent approximation the scattering can 
be treated within ray optics, where rays travel in straight 
lines and only scatter at the line defect, where they will 
either transmit with a probability approximately given 
by Eq. ((!]) or reflect while obeying the law of specular 
reflection. 

The filter can be used to create a valley-polarized beam 
of electrons or holes. By probing the current passing 
through the line defect at a particular angle, one could 
also measure the valley polarization of the incident quasi- 
particles. Demonstration of these components should sig- 
nificantly accelerate research on graphene valleytronics. 

The authors acknowledge support from the U.S. Office 
of Naval Research, directly and through the U.S. Naval 
Research Laboratory. D.G. thanks C.W.J. Beenakker for 
helpful comments. 



[1] K. I. Bolotin, K. J. Sikes, Z. Jiang, G. Fudenberg, 

J. Hone, P. Kim, and H. L. Stormer, Solid State Com- 

mun. 146, 351 (2008). 
[2] A. A. Balandin, S. Ghosh, W. Bao, I. Calizo, D. Tewelde- 

brhan, F. Miao, and C. N. Lau, Nano Lett. 8, 902 (2008). 
[3] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
[4] P. R. Wallace, Phys. Rev. 71, 622 (1947). 
[5] J. Lahiri, Y. Lin, P. Bozkurt, 1. 1. Oleynik, and M. Batzill, 

Nat. Nanotech. 5, 326 (2010). 
[6] P. Y. Huang, C. S. Ruiz- Vargas, A. M. van der Zande, 

W. S. Whitney, M. P. Levendorf, J. W. Kevek, S. Garg, 



5 



J. S. Aldcn, C. J. Hustcdt, Y. Zhu, J. Park, P. L. McEuen, 

and D. A. MuUcr, Nature 469, 389 (2011). 
[7] S. Malola, H. Hakkinen, and P. Koskinen, Phys. Rev. B 

81, 165447 (2010). 
[8] O. V. Yazyev and S. G. Louie, Phys. Rev. B 81, 195420 

(2010). 

[9] Y. Liu and B. I. Yakobson, Nano Lett. 10, 2178 (2010). 



[10] O. V. Yazyev and S. G. Louie, Nat. Mat. 9, 806 (2010). 
[11] A. Rycerz, ,J. Tworzydlo, and C. W. J. Beenakker, Nat. 

Phys. 3, 172 (2007). 
[12] T. Ando and T. Nakanishi, J. Phys. Soc. Jpn. 67, 1704 

(1998). 

[13] J. C. Slonczewski and P. R. Weiss, Phys. Rev. 109, 272 

(1958). 



